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In parameter space of nonlinear dynamical systems, windows of periodic states are aligned fol- 
lowing routes of period-adding configuring periodic window sequences. In state space of driven 
nonlinear oscillators, we determine the torsion associated with the periodic states and identify re- 
gions of uniform torsion in the window sequences. Moreover, we find that the measured of torsion 
differs by a constant between successive windows in periodic window sequences. We call this phe- 
nomenon as torsion-adding. Finally, combining the torsion and the period adding rules, we deduce 
a general rule to obtain the asymptotic winding number in the accumulation limit of such periodic 
window sequences. 



A conspicuous characteristic in parameter space of dis- 
sipative nonlinear dynamical systems is the appearance 
of periodic states for parameter sets immersed in param- 
eter regions correspondent to chaotic states. In the lit- 
erature, much attention has been devoted to establish 
connections between these periodic states. For exam- 
ple, a successive constant increment on the period of os- 
cillation of such states (period-adding phenomenon) [l| 
have been experimental and numerically observed in sev- 
eral real- world systems such as neuronal activities [l, Q , 
electronic circuits bubble formation 0, semiconduc- 
tor device [6[, and chemical reaction fr\. The period- 
adding phenomenon has been also observed for sequences 
of shrimp-shaped periodic windows accumulating in spe- 
cific parameter space regions [8-13|. Once nonlinear dy- 
namical systems can exhibit many different kinds of mo- 
tion, knowing adding rules, such as the period-adding 
and further information about the accumulating parame- 
ter regions, is very advantageous, specially, for predicting 
periodic states for different parameter sets in real-world 
applications. 

Furthermore, besides the intrinsic period of oscilla- 
tions, in dissipative systems, periodic states have other 
interesting convergence properties. For instance, for 
driven nonlinear oscillators, the torsion number n is de- 
fined as number of twists that local flow perform around 
a given periodic solution during a dynamical period m, 
and the winding number defined as w = n/m 
[ist . However, besides the existence of such convergence 
properties, additional connecting rules between periodic 
states and accumulating regions characteristic have not 
yet been discovered. 

Our aim here is to investigate the convergence charac- 
teristics, namely, the torsion and winding number of pe- 
riodic states within complex periodic windows, in period- 
adding sequences in the parameter space of driven non- 
linear oscillators. A torsion-adding formulation between 
such periodic states are proposed here. Combining both 
additive sequences properties, the torsion and the period 
adding, we describe a generic periodic window in a se- 



quence in terms of its winding number. The asymptotic 
limit of such description gives a general rule to deter- 
mine the winding number for any accumulation of period- 
adding sequences. 

Generally, the driven nonlinear oscillator is described 
by: 



X H- g{x, x) = h{t)^ 



(1) 



where h{t) = /i(t + T) is a periodic function with angular 
frequency uj = 27r/T. For this equation, the winding 
number is obtained by considering revolutions performed 
by an orbit around of a very close neighbor periodic 
orbit 7 during the interval time At = T (see Fig. [T]). The 
absolute mean value of the revolution angular frequency, 
^(7) = I < > I? is called torsion frequency: 



1^(7) = lim i 



I a{t')dt' 
Jo 



limMi^i^. (2) 
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Thus, considering the T-shift map [l9|, the winding num- 
ber is precisely defined as Q: 



^(7) = 



(3) 



For a more appropriated form, note that the 7 period is 
given by = mT = 27rm/cj, while the torsion frequency 
period is Tq = 27r/Q (see Fig. [1]). Including these pe- 
riods in Eq. (|3]), and defining the torsion number as 
n = T^/Tq^ we obtain the winding number 



^(7) = — • 
m 



(4) 



Figure 1 here 

In Fig. [21 we show a sketch of the region with the low- 
est period inside a complex periodic window. We refer 
to it as main-body window^ frequently found in the two- 
dimensional parameter space of the dynamical system 
given by Eq. [H This main-body window has the same 
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FIG. 1. (Color online) Sketch of a periodic orbit 7 with period 
surrounded by a neighbor orbit 7^, a{t) is the phase of the 
7^ revolution around 7. The stroboscopic Poincare section S 
is defined by the forcing period T. 



adding phenomenon in periodic window sequences. In 
other words, for each increment p, in the window period 
mi, the torsion number of equivalent regions is also in- 
cremented by a constant value r. Therefore, from Eq. 
m the winding number wr. in a generalized main-body 
window region R {R = A, B^C^ or D) of the i-th window 
sequence can be determined in function of the torsion 
number in region of a known window, in particular, in 
function of the first one, ur-^ : 



riR^ + - 1)t 
mi + (i - l)p ' 



(5) 



Thus, the asymptotic winding number limit of any region 
R, 



period and is bounded by saddle-node (sn) and period- 
double (pd) bifurcation curves. On the other hand, we 
identify that the winding number (Eq. [3]) is the same 
in four different regions A, C, and D. The point 
where the two saddle-node bifurcation curves met in re- 
gion B depicts a cusp bifurcation, where two periodic 
orbits are created. Bellow this point, there are limited 
regions where 5, C, and D appear in pairs {BC^ DB^ 
and DC)^ due to the coexistence of the two periodic or- 
bits. The curves Ai and A2 between the colored regions, 
reveal the window skeleton [13, Hll • The flow converges 
monotonically or non monotonically to the periodic orbit 
according the skeleton composition (Ai, A2). Crossing a 
A curve from one side to the other, the flow convergence 
suffers a transition. The direction of this transition char- 
acterizes Ai and A2 and causes a change in the winding 
number 12211. 
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FIG. 2. (Color online) Main-body window in two-dimensional 
parameter space, sn and pd are the saddle-node and period- 
doubling bifurcation curves. The curves Ai and A2 are the 
skeleton of the structure delimiting the regions A (red), B 
(gray), C (yellow), and D (green) are regions delimited by 
the skeleton. 

Figure 2 here 

Then, we consider the winding number concept for se- 
quences of periodic windows. Until now these sequences 
have been described only by period-adding rules, i.e., the 
period of a periodic window in the sequence can be deter- 
mined by adding a constant value p to the period of the 
previous window 23l-[26|. Now, we introduce the torsion- 



shows that the winding number of all regions converge 
to a constant which only depends on the torsion number 
and the period increments, r and p, respectively. We 
denote this limit as Woq. 

Another consequence of the torsion-adding is that the 
skeleton of all windows in a sequence are equivalent. In 
fact, one can show that for Ai curve, in successive win- 
dows i and i + 1, the torsion-adding condition {nR._^^ = 
ur. +r) imply that the torsion number difference between 
regions A and C is the same (n^^ —Tid = ^Ai+i —^d+i)- 
Thus, the Ai curve promotes the same transition in all 
sequence. The same is for A2. 

Now, to verify our results, we present numerical simu- 
lations for a specific driven nonlinear oscillator, namely, 
the Morse oscillator that describes a diatomic molecule, 
immersed in an external electromagnetic field, modeled 
by the following differential equation [27|: 



dx 



') = 2.5cos(cjt), 



(7) 



where the parameter d is the amplitude of the system 
damping and uj is the angular frequency of the external 
forcing. Our analysis are in the two parameter space 
d X CO. 

As we are interested in how a trajectory 7' converges to 
the stable periodic orbits 7, we consider in our numerical 
simulation 7' starting in a position very close to 7 in a 
such way that the linearized flow is enough to describe 7' 
dynamics. We represent a point in 7 as {xl^X2,x^) and 
in y as (^1,^2,^3)- Therefore, the 7' evolution is given 

by 



yi = y2 

y2 = 8e-<{l - 2e-<)yi - dy2, 



(8) 



where we considered ^3 = 0, since ys is any constant 
in the linearized flow, i.e., ^3 = 0. Finally, considering 
the solutions of Eqs. (j8|) in the polar coordinates yi = 



3 



rcos(a) and y2 = rsin(a), we determine the angular 
frequency, 



winding number is given by 



-vl 



(2/12/2 - 2/22/1 )■ 



(9) 



Thus, by Eqs. (|2]) and (|3]), we can determine numeri- 
cally the winding number inside a period-m main-body 
window and then calculate the associated torsion number 
with Eq. dlj, where m is given hy m = T^/T. 

We recall that the orbit 7 is stable if the Floquet multi- 
pliers jiii associated with the mT-shift map (or any multi- 
ple of mT) are < 1. For jUi G M, if /ii are positive, the 
orbit 7' converge to 7 with orientation preserving (mono- 
tonically) and if fii are negative, 7' converge to 7 with 
orientation reversing (non monotonically) [22]. Since the 
multipliers are < 1 from the saddle- node bifurcation 
(/i = 1), in region A, to the period-doubling bifurcation 
(/i = —1) in regions C and the phase difference by 
the interval time At = mT computed between regions A 
and C (or D) is tt or, i.e., \nA — nc{D) \ = 1/2, as can be 
seen in Fig. [3l 

For the Morse oscillator, we obtain two-dimensional 
parameter spaces [shown in Figs. [3]and|4] by computing 
the torsion and the winding numbers for each d x uj pa- 
rameters of a two-dimensional mesh of 500 x 500 equally 
spaced. We assign different colors to designate wind- 
ing number values of periodic attractors, while the white 
color represents the parameters correspondent to chaotic 
attractors. With this procedure we identify the uniform 
winding number regions, A, 5, C, and D indicated in 
Fig. [2j We also identify inside the three periodic win- 
dows, in Fig. [31 the skeleton separating areas with uni- 
form torsion numbers. To be more precise, lets define 
the left and right side of the curves Ai and A2 traversing 
the curve from the bottom to the top (see Fig. [2]). The 
torsion number is always increased or decreased by 1/2 
when we cross the curve Ai (A2) from the right (left) to 
the left (right). The superscripts + and — , in Ai and A2, 
indicate if the torsion number increases (+) or decreases 
(— ) according to the defined orientation. All the three 
possible skeleton are shown in Fig. [3l 

Figure 3 here. 

Note that the main-body windows, shown in Fig. [3l 
presents three different winding number values. In Fig. 
[3ja), the torsion number decreases by 1/2 from A to C 
(Ai = A^, according our convention) and increases by 
1/2 from C to B (A2 = A2"). Thus, the torsion number 
difference between regions A and B is zero. The same is 
verified for the route A ^ D ^ B^ then ua = tib- Since 
the parameter sets in this window correspond to orbits 
with the same period m, according Eq. (|5]), the winding 
number in regions A and B are the same. Similarly, we 
conclude that in Fig. ^h) and (c), where nc = ^d, 
the winding number of regions C and D are the same. 
Therefore, for the i-th complex window (see Fig. H]), the 



nAi + (i-l)T + (fe-/)/2 
mi + (i - l)p 



(10) 



where k^l G {0,1,2} are, respectively, the number of 
times that the curves type A+ and A~ are crossed from 
region A to any region R crossing Ai and A2 just one time. 
Equation (fTQ|) describes the internal regions in these pe- 
riodic window sequence and obey the convergence limit 
established in Eq. (|6]). Moreover, Eq. (fTQ|) states that 
if one periodic windows presents wb = t/p, then, the 
winding number of region B in any periodic window in 
the sequence is WBi = r/p. 

To illustrate the validity of our results, we display in 
Fig. m^a) sequences of periodic windows where all win- 
dows have their internal regions separated by A^ and 
A J. For this sequence, according to Eq. (fTQ|) with k = 1 
and / = 1, all periodic windows present different winding 
number values in their central region B. Additionally, we 
show in Fig. [H^b) the winding number calculated along 
a line passing through the region B of all windows that 
compose the sequence. It is clear that the winding num- 
bers converge to Woo = 3.0. The torsion number and the 
period increment can also be determined, in Fig. H^a), 
by r = n^.+i - tib, = 12 and p = m^+i - = 4, re- 
spectively. Thus, the winding number convergence is in 
agreement with Eq. (|6]). 

Figure 4 here. 

In Fig. Sl^c), we show another sequence of periodic 
windows internally separated by and A J. For this 
sequence of windows, we measure r = 5, p = 3, and 
the winding number wb = 5/3. Thus, as predicted by 
Eqs. (p!Q|) with k = 2 and / = 0, we measure WBi =5/3. 
From Fig. [41(d) we also verify the limit Woo = r/p = 
5/3. Figure [4] shows that the regions where the sequences 
are accumulating have the same correspondent winding 
number w^q. 

In conclusion, we report the existence of torsion-adding 
phenomenon in periodic window sequences of the driven 
nonlinear oscillators. Additionally, we formulate a gen- 
eral rule [Eq. ([TQ|) ] to obtain the winding number of any 
window belonging to the sequence. From this general rule 
we obtain the winding number asymptotic limit for any 
sequence {woo = ^/p), such ratio seems to be an univer- 
sal property of dynamical systems, once it requires only 
the existence of the period and torsion adding phenom- 
ena. Moreover, since there are no general theory ensuring 
that period-adding sequences are composed of an infinite 
number of windows, this limit is a theoretical evidence of 
the existence of infinite windows in the sequences. 

Furthermore, we present numerical simulations for the 
Morse oscillator where the reported torsion-adding phe- 
nomenon and winding number asymptotic limit are veri- 
fied. We also performed numerical analysis for other non- 
linear oscillators described by Eq. ([1]), and verified that 
all results are in complete agreement with our theory. 
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FIG. 3. (Color online) The white region indicates chaos while the colored regions indicate complex periodic windows according 
to the winding numbers of a grid of the system parameters. The numbers noted on the different winding number regions, A,B, 
C, and D, are the correspondent torsion numbers n. The A symbols indicate the window skeleton, and the A superscript + or 
— indicate increases or decreases in the torsion numbers, (a) ua = ns = 5.0, nc = 5.5, and no = 4.5. (h) nc = no = 14.5, 
ua — 14.0, and ub — 15.0. (c) nc — tid — 15.5, ua — 16.0, and ub — 15.0. 




FIG. 4. (Color online) (a) Periodic window sequences obtained by computing the winding numbers. The pairs (nB,ms) noted 
on different periodic windows are the torsion number and the period of region B. The increments on the torsion number and 
periods are, respectively, r = 12.0 and p — 4.0. (b) Winding number, as a function of parameter cj, calculated along a line 
crossing the B regions of the window sequence of (a). Here, the asymptotic winding number for this sequence is Woo = 3.0, 
equal to its predicted value r/p = 3.0. (c) The same as (a) in another region of the parameter space (in this case r = 5.0 and 
p = 3.0). (d) The same as in (b) (in this case Woo = r/p = 1.6666). The colors are the same meaning of Fig. (3] 
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